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Abstract  -  The  operation  of  a  switched  reluctance  machine 
with  eccentric  rotor  positions  creates  asymmetrical  airgap  flux 
distributions  and  results  in  unbalanced  magnetic  pull.  This  paper 
comprehensively  investigates  the  static  and  dynamic  radial  force 
characteristics  of  an  18/12  switched  reluctance  machine  with  un¬ 
balanced  operation.  A  powerful  magnetic  equivalent  circuit 
(MEC)  modelling  approach  is  developed  and  allow  fast  and  accu¬ 
rate  force  predictions  with  arbitrary  excitation  and  position.  The 
model  calculates  the  instantaneous  radial  and  tangential  forces  on 
the  individual  poles  as  wellasthe  net  forces  on  the  rotor.  The  meth¬ 
od  is  validated  with  finite  element  analysis  (FEA)  and  against  ex¬ 
perimental  results  from  a  12/8  machine.  Static  radial  force  charac¬ 
teristics  are  first  developed.  The  MEC  is  then  coupled  with  a  time 
stepping  algorithm  to  calculate  the  simultaneously  varying  electri¬ 
cal,  mechanical,  and  magnetic  variables  for  dynamic  operation  of 
the  complete  SR  motor  drive.  Dynamic  radial  force  behavior  is  ex¬ 
amined,  where  parallel  connected  windings  are  found  to  naturally 
and  significantly  reduce  the  unbalanced  magnetic  pull.  Variations 
of  radial  force  with  torque-speed  operating  point  are  also  ex- 
amtned. 


I.  Introduction 

Switched  Reluctance  Machines  (SRM’s)  have  been  the 
subject  of  much  active  research.  Presently  they  are  finding 

their  way  into  various  commercial,  industrial,  and  military  ap¬ 
plications  [1].  Perhaps  one  of  the  most  established  applica¬ 
tions  is  that  of  a  starter/generator  for  an  advanced  military  mr- 
engine.  Specific  reejuirements  of  fault  toleranc^e,  high 
temperature,  and  high  speed  make  the  SR  machine  a  strong 
candidate  for  application  to  aircraft  power  generation  system 
[2)(3][4].  Additional  advantages  associated  with  high  speed 
operation,  low  friction,  and  active  positioning  are  realized 
when  the  machine  is  integrated  with  a  magnetic  bearing  sys¬ 
tem  and  located  direcnly  on  the  aircraft  engine  shaft. 

The  potential  advantages  of  SR  technology  present  re¬ 
searchers  and  engineers  with  some  exciting  technical  chal¬ 
lenges  in  the  areas  of  machine  design,  controls,  power  elec¬ 
tronic  drive  circuitry,  analysis  and  system  integration  to  name 
a  few.  These  areas  have  been  the  subject  of  much  research  and 
development.  While  many  of  these  areas  are  the  subject  for 
ongoing  researchl  little  has  been  done  to  address  the  implica¬ 
tions  of  integration  with  a  magnetic  bearing  system.  The  mag¬ 
netic  bearing  system  presents  an  additional  5  degrees  of  free¬ 
dom  for  the  mechanical  rotor  motion,  thus  requiring  eccentric 
operation  of  the  SR  machine.  The  eccentric  positioning  of  the 
SR  rotor  creates  unbalanced  electromechanical  radial  forces 
due  to  asymmetrical  magnetic  pull.  A  thorough  understand¬ 
ing  of  unbalanced  radial  forces  and  detailed  models  that  can 


predict  them  for  various  modes  of  operation  with  various  ma¬ 
chine  geometry  is  required  and  could  be  most  useful  for  study¬ 
ing  other  conventional  SR  drive  issues,  such  as  the  impact  of 
mechanical  and  manufacturing  imperfections  on  the  load  and 
hence  life  of  conventional  bearings,  and  acoustic  noise  issues. 

The  principle  objective  of  this  research  is  to  characterize 
the  radial  force  properties  of  the  SR  machine  associated  with 
eccentric  rotor  positiom'ng.  Radial  forces  must  be  predicted 
and  understood  for  establishing  the  bearing  performance  re¬ 
quirements  as  well  as  providing  a  means  for  compensation  of 
unwanted  disturbances  to  both  radial  and  angular  dynamics. 
To  this  aim,  a  versatile  magnetic  equivalent  circuit  (MEC) 
model  is  developed  for  the  purpose  of  characterizing  the  ma¬ 
chine  behavior  under  a  wide  range  of  operating  modes.  The 
model  is  required  to  enable  fully  coupled  nonlinear  magnetics, 
arbitrary  excitation,  and  to  be  capable  of  predicting  the  elec¬ 
tromagnetic  field  and  electromechanical  force  production 
with  unbalanced  mechanical  operation.  Addidonal  capability 
is  required  to  couple  with  the  temporal  differential  equations 
of  electrical  excitation  and  mechanical  motion  to  provide  a 
fully  dynamic  simulation  tool.  A  series  of  static  force  inves¬ 
tigations  are  performed  to  characterize  the  radial  force  prop¬ 
erties  as  a  function  of  instantaneous  position  and  excitation. 
Particular  to  this  study,  is  the  motivation  to  examine  the  dy¬ 
namic  radial  force  properties  for  various  modes  of  operation. 
The  results  can  then  be  extended  to  quantify  the  dynamic  stiff¬ 
ness  contributed  to  the  entire  rotating  assembly  by  the  unbal¬ 
anced  SR  machine.  The  models  are  developed  to  flexibly 
cope  with  various  SR  machine  geometries  such  as  6/4, 12/8, 
and  18/12.  The  results  of  this  effort  specifically  address  an 
18/12  250kw  machine  and  provide  a  comprehensive  character¬ 
ization  of  the  imbalanced  radial  forces  for  various  mechanical 
and  electrical  operating  modes.  The  parameters  for  the  ma¬ 
chine  under  investigation  are  given  in  Thble  I. 

TABLE  I 

Key  Parameters  of  the  18/12  Switched  Reluctance  Machine 


Parameter 

Vilue 

Stator  outside  diameter 

33.^icm  (13.36’^) 

Rotor  outside  diameter 

25.6cm  (lO.n 

Stack  length 

11.4cm  (4.51 

Airgap  length 

0.762mm  (0.030") 

Maximum  expected  eccentricity 

33%  of  airgap 

Lamination  material 

HyperooSO 

Rated  generating  power 

250  kW 

Maximum  operating  speed 

18,000  rpm 

Rated  DC  bus  voltage 

270  V 

Rated  peak  pole  current 

250  A 
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II.  Magnetic  Equivalent  Circuit  Modelling 
A.  Circuit  Structure 

A  powerful  magnetic  equivalent  circuit  (MEC)  based 
model  is  developed  to  provide  a  versatile  means  for  investigat¬ 
ing  the  unbalanced  SR  machines  of  this  study.  The  need  for  an 
efficient,  accurate  and  static  as  well  as  dynamic  model  thatc.^ 
cope  with  the  complexity  of  an  entire  machine  with  meclm^ 
cal  or  electrical  a^metry  motivates  the  selection  of  MEL 
analysis.  Magnetic  Equivalent  circuit  modeling  affords  a  good 
compromise  between  the  accuracy  of  finite  element  and  the 
computational  speed  of  analytical  models.  The  MEC  nemork 
is  designed  from  the  machine  geometry  and  from  a  fore¬ 
knowledge”  of  the  flux  paths.  The  circuit  consists  of  a  network 
of  passive  permeance  elements  and  active  M^  source  ele¬ 
ments.  The  passive  permeance  elements  consist  of  both  satu¬ 
rating  iron  and  linear  air  elements  both  of  fixed  and  variable 
geometries.  The  iron  magnetic  characteristics  are  govei^d 
by  the  inherent  BH  properties  of  the  selected  matenal.  The 
result  is  a  non-linear  magnetic  circuit  including  mutual  wu- 
pling  with  independent  and  arbitrary  excitation  and  position¬ 
ing.  The  model  canbe  embedded  in  a  time  stepping  algorithm 
for  simulation  of  dynamic  behavior  coupled  at  the  mechamcal 
and  electrical  terminals.  Fig.  1  illustrates  a  simplified  MEC  for 
a  6/4  SR  machine.  For  convenience  of  illustration,  only  core, 

tooth  and  gap  flux  paths  are  shown  on  this  diagram.  Note  that 

there’is  individual  and  independent  pole  excitation  and  that 
there  is  no  symmetry  assumed  as  with  most  SR  models. 


Fig.  1  Simplified  MEC  for  6/4  SRM 
B.  Computational  Engine 

The  computational  implementation  forthe  model  is  devel¬ 
oped  to  provide  a  maximum  level  of  versatility.  Tb  this  aim,  a 
generic  formulation  is  developed  to  automatically  deal  with 
different  machine  configurations  and  operating  modes,  in¬ 
cluding  geometry,  stator/rotor  pole  number,  windings,  mag¬ 
netic  material,  and  excitation.  The  computational  structure  is 
comprised  of  several  key  parts  based  on  the  MEC  form  of  Fig. 
1  A  detailed  definition  of  the  machine  parameters  and  input 
structure  initializes  the  analysis.  Next  a  set  of  matrix  equations 
are  automatically  constructed  whose  coefficients  and  inputs 
are  calculated  from  the  instantaneous  excitation  and  position 
information.  Anodal  analyas  approach  is  selected  for  its  gen¬ 
erality,  rigor  and  ability  to  handle  non-planar  dremts,  unlike 


a  mesh  based  dreuit  approach.  A  numerical  solution  algo¬ 
rithm,  forming  the  heart  of  the  computational  model,  itera¬ 
tively  calculates  the  instantaneous  non-linear  magnetostatic 
solution.  A  post  processing  algorithm  calculates  the  desired 
output  variables  including  individual  radial  and  tangential 
forces  for  each  pole,  net  force  on  the  rotor,  flux  and  MMF  dis¬ 
tributions,  flux  linkages,  etc..  The  Maxwell  stress  tensor  com¬ 
bined  with  an  appropriately  selected  integration  path  easily 
enable  the  force  calculations  for  each  pole.  Coordinate  trans¬ 
formations  are  utilized  to  vectorially  sum  the  radial  forces  to 
determine  the  net  radial  force  acting  on  the  rotor  body. 

III.  Finite  Element  Modelling 

An  automated  2-D  finite-element  analysis  (FEA)  pack¬ 
age  for  switched  reluctance  machines  was  prepared  by  a  finite 
element  software  company.  The  package  is  enormously  flex¬ 
ible  in  evaluating  magnetostatic  situations.  Rotor  and  stator 
grids  are  automatically  created  from  a  data  file  of  the  ma¬ 
chine’s  major  dimensions.  The  angular  coordinate  of  the  rotor 
and  the  horizontal  &  vertical  off-center  deflections  are  used 
to  position  the  rotor  grid  with  respect  to  the  stator  gnd.  The 
current  in  each  stator  coil  is  individually  specified,  and  this, 

combined  with  rotorpositioningflexibility,  allowsamagnetos- 

tatic  evaluation  of  any  possible  situation  of  an  SR  machine. 
Automatic  post-processing  calculates  the  electromagnetic 
torque  and  the  horizontal  &  vertical  electromagnetic  forces 
from  a  magnetostatic  solution  of  the  machine.  The  torque  and 
force  for  each  individual  rotor  pole,  as  well  as  a  summed  total 
rotor  torque  and  force  are  calculated.  Fig.  2  illustrates  the  flux 
plot  of  an  18/12  SRM  with  one  phase  excited  for  an  offeet  rotor 
position.  Note  the  greater  flux  density  on  the  right  side  of  the 
machine  where  the  airgap  is  reduced,  resulting  in  an  unbal¬ 
anced  magnetic  pull. 
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Fig.  2  FEA  Flux  plot  for  an  1 8/12  SRM  with  a  33%  gap  rotor  offset 
position. 

IV  Flux  Linkage  Characteristics 

The  flxix  linkage  characteristics  of  the  SR  machine  form 
the  heart  of  virtually  all  analytical  SR  models.  The  flux  linkage 
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characteristics  indicate  much  of  the  principle  behaviors  of  th 
machine,  namely  inductance  and  torque  production.  U^ike 
the  analytical  models,  the  MEC  like  FEA  has  the  capability  to 
generate  the  flux  linkage  curves  of  the  machine.  It  is 
fitting,  as  a  first  step  in  verification  and  calibration  of  the 
model  to  calibrate  the  model  with  based  on  the  flux  linkage 
curves’of  the  machine.  The  finite  element  model  is  med  to 
generate  a  reference  set  of  curves  for  comparison.  Fig.  3  il  us- 
trates  the  flux  linkage  characteristics  of  the  machine  calcu¬ 
lated  with  rotor  position  as  a  parameter  for  balanced  opera¬ 
tion  Note  that  balanced  operation  here  refers  to  symmetrica 
electrical  excitation  of  the  poles  and  to  a  centered  mechamc^ 
rotor  position,  i.e.  no  eccentricity.  It  is  obseived  that  the 
predictions  agree  very  well  with  those  obtained  from  F^. 
This  illustrates  the  capability  of  the  MEC  model.  Most  SRM 
modeling  techniques  require  the  flux  linkage  charactcnstics  a 
an  input,  whereas  the  MEC  model  has  inherent  capability  to 
predict  the  characteristics.  It  should  also  be  noted  that  in  geri- 
eral  the  flux  linkage  characteristics  are  a  function  of  four  van- 
ables  rather  than  two  for  the  balanced  machine.  The  eccen¬ 
tricity  contributes  two  additional  variables  to  the  usual 
variables  of  current  and  angular  position.  Equation  (1)  mdi- 
cates  the  dependency. 


A  =  A(i,0,A:,y) 

The  MEC  model  must  be  capable  of  dealing  with  the  add¬ 
ed  complexity  to  accurately  model  the  flux  distributions  and 
hence  the  radial  force.  In  addition  to  the  eccentricity,  the 
poles  are  able  to  be  independently  excited  as  an  added  com¬ 
plexity.  Typically  SR  analyses  assume  symmetrical  exatation 
as  well  as  geometric  symraetiy  and  are  incapable  of  dealing 
with  unbalanced  and  unsymmetrical  operation.  The  capabili¬ 
ty  of  the  MEC  model  to  deal  with  these  complexities  also  make 
it  suitable  for  the  investigation  of  faulted  operation,  a  subject 
of  much  interest  for  the  SR  machine. 


Fig.  3  Flux  Linkage  Characteristics,  18/12  SRM  -  Comparison  of  MEC 
&FEA 


V.  Static  Radial  Force  Characteristics 

A.  Radial  Force  Computation 

Before  presenting  the  predicted  radial  force  characteris¬ 
tics,  some  explanation  of  the  radial  force  calculaUon  is  in  or¬ 
der!  The  calculation  of  radial  forces  are  based  on  the  Maxwell 
Stress  Tensor  in  tlic  developments  of  this  work.  Others  [5] 
have  presented  energy  based  methods  involving  the  differen¬ 
tial  variations  in  the  geometric  dependent  permeance  ele¬ 
ments,  which  can  be  shown  equivalent.  Since  the  MEC  by  na¬ 
ture  assumes  flux  paths  apriori  to  the  solution,  the  assumed 
flux  lines  are  well  defined  and  very  suitable  for  application  of 
the  Maxwell  stress  tensor.  Together  with  the  appropriate 
selection  of  integration  paths  around  the  individual  pole  faces, 
the  Maxwell  stress  tensor  provides  straight  toward  calcula¬ 
tion  of  the  radial  forces  once  the  magnetostatic  solution  is  ob¬ 
tained.  Recall  the  Maxwell  stress  tensor  as  given  in  (2),  where 
the  resultant  force  is  independent  of  the  integration  path. 


B\)dA 


(2) 
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To  illustrate  consider  the  alignment  of  two  poles  as  indi¬ 
cated  in  simplistic  cartoon  representation  of  Fig.  4.  The 
approximations  of  the  flux  paths  being  always  radial  across  the 
airgap  ease  the  selection  and  calculation  of  the  stress  contour. 
The  contour  is  selected  such  that  the  flux  lines  are  always  ei¬ 
ther  orthogonal  or  parallel  to  the  contour .  Under  these  condi¬ 
tions  the  Maxwell  stress  equation  simplifies  to  yield  only  force 
normal  to  the  contour  as  follows: 


Fig.  4  Application  of  the  Maxwell  Stress  Tensor  in  the  calculation  of 
radial  and  tangential  forces. 


Orthogonal  Field  Lines:  -* 

Parallel  Field  Lines:  =  ^|  | 

Applying  these  equations  to  the  contour  illustrated  in  Fig. 
4  residts  in  the  following  radial  and  tangential  forces  acting  on 
the  rotor  tooth: 


(4) 


(5) 
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where  Bn,  Bf2,  and  Bm  are  the  first  and  second  airgap  fringing 
fields  and  the  main  airgap  field,  respectively.  Thus  the  radial 


force  is  equivalent  to  integrating  surface  area  of 

dW 

the  airgap  and  the  tangential  force  is  equivalent  to  where 

ITis  the  stored  energy  in  the  overlapping  region  of  the  airgap. 
Equation  (7)  and  the  rotor  radius  are  directly  used  to  calculate 
the  torque  as: 

T -- 

Lastly,  all  of  the  force  contributions  of  each  active  pole 
pair  are  combined  together  to  determine  the  net  radial  force 
acting  on  the  rotor  body.  This  entails  accounting  for  the  onen- 
tation  and  magnitude  of  each  force.  A  coordinate  transforma¬ 
tion  local  to  each  pole  pair  easily  maps  the  local  forces  to  the 
main  x,y  stator  reference  system.  The  forces  and  torques  are 
then  vector  summed.  The  coordinate  transformation  for  one 


sudi  pole  pair  is  given  as: 


In  summary,  the  appplication  of  the  Maxwell  stress  tensor 
to  the  well  defined  magnetic  circuit  flux  paths  is  appealing. 


simple  and  intuitive. 

B.  Force -vs- Pole  Current 

A  fundamental  characteristic  of  the  radial  forces  in  the  ma- 
chine  as  a  function  of  current  and  eccentricity  is  the  force  ver¬ 
sus  pole  current  with  eccentricity  as  a  parameter.  Fig.  5  illus¬ 
trates  the  force— vs-pole  current  for  the  18/12  SR  machine  for 
two  different  eccentricities  with  the  machine  in  the  fully 
aligned  positioa  The  plot  also  compares  the  predictions  of 
the  lumped  circuit  model  vrith  that  of  finite  element.  It  is  ob¬ 
served  that  the  MEC  predicts  well  with  the  FEA.  Charactens- 
tically,  the  force  is  observed  to  rise  with  the  square  of  current  in 
the  linear  region  as  expected,  reach  a  peak  and  then  fall  off 
dramatically  wth  magnetic  saturation.  The  effect  of  eccen¬ 
tricity  is  also  observed  from  the  plot  where  an  increased  rotor 
displacement  produces  a  larger  unbalanced  radial  foree.  In 
the  linear  region,  flux  builds  linearly  with  current  and  hence 
the  force  goes  as  the  square.  As  the  knee  is  reached  and  the 
machine  begins  to  saturate,  the  force  peaks.  As  saturation  sets 
in  the  flux  tends  to  equalize  around  the  machine  and  the  un¬ 
balance  dissipates.  In  hard  saturation,  the  radial  force  drops 


off  dramatically.  Another  view  is  that  in  the  extreme,  the  satu¬ 
ration  effectively  widens  the  airgap  making  the  eccentricity 
much  effective  smaller  and  the  reluctance  is  balanced  around 
the  machine.  As  expected  the  largest  discrepancies  between 
the  two  predictions  is  in  the  saturated  region.  Saturation  is  the 
most  difficult  area  to  accurately  match  with  MEC  as  the  Dux 
paths  become  much  less  defined. 

18/12  SRM  100%  Alignment 


Pole  Current  (amp) 

Fig.5  Radial  force  -vs-  pole  current,  18/12  SRM,  25%  &  33%  gap 
offset  in  rotor  position. 

It  is  important  to  recognize  that  the  net  radial  force  is  the 
unbalance  of  large  radial  forces  on  a  per  pole  basis.  Fig.  6  illus- 
trates  the  radial  force  contributions  of  the  two  diametrically 
opposite  poles  in  line  with  the  eccentricity.  The  forre  charac¬ 
teristics  of  Fig.  5  are  more  apparent  when  the  individual  pole 
forces  are  considered. 

18/12  SRM  100%  Alignment 


Fig.  6  Radial  fpree  per  pole  -vs-  pole  current,  18/12  SRM,  25%  & 
33%  gap  offset  in  rotor  position. 


It  should  be  noted  that  each  point  on  these  curvw  repre¬ 
sents  a  potential  instantaneous  point  under  dynamic  opera¬ 
tion.  Also  the  peak  force  here  of 800  lbs  represents  the  largest 
instantaneous  force  obtained  for  operation  limited  to  this  ec- 
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centricity  and  peak  current  level.  It’s  also  interesting  to  keep 
in  mind  that  the  worst  case  radial  force  does  not  occur  at  maM- 
mum  excitation.  These  results  will  be  revisited  again  in  the  dis¬ 
cussion  of  dynamic  operation. 

C.  Force  -vs-  Angular  Position 

Another  fundamental  characteristic  is  that  of  the  radial 
force-vs-pole  position  vnth  eccentricity  and  current  as  pa¬ 
rameters.  Fig.  7  illustrates  the  radial  force  for  a  33  %  eccentric¬ 
ity  as  a  function  of  rotor  angular  position  with  current  as  a  pa¬ 
rameter.  In  this  case  one  set  of  poles  or  one  phase  is  exated 
with  a  fixed  DC  current  and  then  the  rotor  is  rotated  from  the 
unaligned  position,  through  alignment  and  into  the  next  un¬ 
aligned  position.  It  is  observed  that  the  peak  forces  occur  at 
100%  alignment  as  expected  from  Maxwell  stress  and  maxi¬ 
mum  flux  density  in  the  gap.  The  affect  of  varying  current  is 
also  observed  in  the  variation  of  the  radial  force  levels.  The 
maximum  force  occurs  at  the  peak  force -vs-current  levels 
and  100%  alignment  as  observed  on  the  previous  plots.  To¬ 
gether  this  plot  and  the  previous  plot  indicate  that  the  maxi¬ 
mum  instantaneous  radial  force  will  occur  at  the  aligned  posi¬ 
tion  with  a  pole  current  of  90  amps.  This  is  to  be  compared 
with  the  rated  peak  pole  current  pf  approximately  250  amps. 


Fig.  7  Radial  forcc-vs— pole  position,  18/12  SRM,  33%  gap  offset  in 
rotor  position,  FEA  results. 


D.  Force  -vs~  Eccentricity 

Lastly  to  characterize  the  static  radial  force  properties  of 
the  machine  is  the  variation  of  radial  force  with  ecxentricity.  In 
this  case  the  current  and  the  the  alignment  position  are  held  at 
the  peak  producing  values,  90  amps  and  100%  alignment,  re¬ 
spectively.  As  expected,  the  unbalanced  radial  force  grows  al¬ 
most  linearly  with  increasing  eccentricity.  At  larger  displace¬ 
ment  the  rate  of  increase  in  force  drops  off.  This  drop  off  is 
most  likely  due  to  the  onset  of  saturation  in  the  magnetics  with 
the  decreased  airgap. 


Eccentricity  (%  airgap) 


Fig.  8  Radial  force -vs -eccentricity,  18/12  SRM,  90  amps  and  100% 
alignment. 


VI.  Dynamic  Radial  Force  Predictions 

A.  Dynamic  Modeling 

So  far  only  static  characterizations  of  the  SR  radial  forces 
have  been  discussed  for  spedfic  operating  points  of  position 
and  excitation.  It  is  however  of  great  interest  to  understand 
the  radial  force  behavior  that  occurs  as  the  machine  is  oper¬ 
ated  dynamically.  This  includes  commutation,  regulation  and 
general  operation  over  the  torque-speed  operating  range. 
While  the  finite  element  program  developed  is  most  accurate 
for  static  calculations,  it  does  not  perform  dynamic  analyws. 
In  order  to  perform  dynamic  time  domain  simulations,  a  time 
stepping  finite  element  program  would  be  required.  Although 
the  results  would  be  very  accurate,  the  already  long  static  com¬ 
putation  time  would  become  unacceptable  for  the  present 
computing  capabilities  available,  espedally  for  the  numerous 
and  varied  cases  considered  for  an  integrated  magnetic-bear¬ 
ing/starter— generator  system.  The  lumped  magnetic  equiva¬ 
lent  circuit  model  (MEC  model)  is  much  faster  with  minimal 
compromise  in  accuracy  and  can  be  coupled  to  the  temporal 
electrical  and  mechanical  differential  equations.  This  enables 
time  domain  simulation  with  unbalanced  mechanic^  and/or 
electrical  conditions  and  dynamic  prediction  of  radial  forces. 
To  this  end,  algebraic  equations  of  the  static  model  are 
coupled  to  the  temporal  differential  equations  of  the  electrica 
and  mechanical  terminal  variables  and  embedded  into  a 
time-stepping  algorithm  to  form  a  complete  dynamic  system 
simulation.  Fig.  9  illustrates  the  MECmodel  embedded  into  a 
dynamic  simul  ation  model  including  the  SR  commutation  and 
torque  control  algorithms,  and  the  mechanical  dynamics.  This 
forms  the  most  powerful  aspect  of  the  MEC  model,  its  ability 
to  combine  computational  speed  with  modelling  accuracy  and 
predict  the  simultaneously  varying  dynamic  behavior  of  the 
electrical,  mechanical,  and  magnetic  variables  coupled  to  the 
entire  switched  reluctance  drive  system.  To  the  authors’ 
knowledge,  no  other  model  of  this  capability  has  been  devel¬ 
oped  for  the  SR  machine. 
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Fig.  9  Dynamic  simulation  model  of  the  SR  system. 

B.  Current  Source  Excitation 

The  first  dynamic  case  to  be  considered  is  that  of  current 
source  excitation.  In  this  case  an  ideal  current  source  is  as¬ 
sumed  and  the  poles  of  each  phase  are  forced  to  have  equal 
currents.  Since  the  magnetic  circuits  will  not  be  symmetrical 
due  to  the  eccentricity,  a  flux  distortion  and  hence  unbalanced 
radial  forces  are  expected  around  the  circumference  of  the 
machine.  Current  excitation  is  approximated  by  low  speed  op¬ 
eration  of  the  machine  where  the  current  regulation  is  very 
strong.  It  is  also  approximated  by  the  case  where  the  poles  of 
each  phase  are  coimected  in  series,  thereby  forcing  the  cur¬ 
rents  to  be  the  same  in  each  pole.  Current  source  excitatioii  is 
a  worst  case  scenario  as  in  general  voltage  constraints  require 
the  poles  to  be  coimected  in  some  parallel  fashion  which  ex¬ 
hibits  inherent  flux  balancing.  To  further  validate  the  force 
calculations,  the  current  waveforms  are  analyzed  statically  on 
a  point— by — point  basis  with  the  FEA  model.  A  current  wave¬ 
form  as  a  function  of  position  or  time  obtained  from  a  simula¬ 
tion  or  real  data  is  forced  into  each  winding  equally  for  a 
point-by-point  calculation  with  FEA.  This  gives  somewhat 
of  a  quasi-dynamic  FEA  solution.  Fig.  10  illustrates  the 
forces  as  calculated  with  FEA  and  the  MEC model.  This  is  the 
first  calculation  of  this  paper  shown  to  account  for  the  forces 
contributed  spatialfy  for  all  of  the  poles  and  phases  as  they  are 
sequentially  excited.  The  plot  illustrates  very  good  agreement 
between  the  FEA  and  MEC  predictions.  Observe  that  the 
force  pulsates  at  a  rate  equal  to  the  number  of  rotor  poles 
times  the  number  of  phases  times  the  synchronous  speed 
(3xl2xRPM  =  7.8kHz  in  this  case,  and  10.8kHz  at  the  maxi¬ 
mum  speed  of  ISkRPM).  The  radial  force  is  also  observed  to 
have  an  average  or  DC  value.  The  force  is  virtually  all  in  the 
direction  of  the  eccentricity.  As  shown  they-directed  force  is 
very  low  in  amplitude  with  virtually  no  DC  value.  The  average 
radial  force  in  the  x-direction  Is  303  lbs  for  the  FEA  data  and 
307  lbs  for  the  MEC  data.  It  should  be  noted  that  each  inst^- 
taneous  point  on  the  plot  in  Fig.  10  corresponds  to  a  condition 
represented  in  the  series  of  static  characteristics  previously 
presented. 

C.  Voltage  Source  Excitation 

A  more  realistic  case  to  be  investigated  is  that  of  voltage 
mode  excitation  with  parallel  connected  poles.  This  is  the  true 


18/12  SRM,  13000  RPM,  250kW 


Fig.  1 0  Radial  force  with  dynamic  current  source  excitation,  1 8/1 2  SR 
machine  with  33%  gap  offset. 

physical  means  by  which  the  machine  will  be  excited.  Even 
with  current  regulators,  the  applied  forcing  function  to  the  ma¬ 
chine  terminals  is  still  voltage.  This  dynamic  analysis  can  only 
be  performed  with  the  MEC  model.  To  illustrate  the  influence 
of  parallel  windings  and  the  capability  of  the  MEC  m^el,  the 
same  simulation  is  repeated  for  voltage  mode  excitation.  For 
voltage  excitation,  the  voltage  and  hence  approximately  the 
fluxes  are  equal  and  the  currents  distribute  accor^ng  to  induc¬ 
tance  or  relative  airgap.  This  is  a  self  compensating  afiect  that 
tends  to  reduce  the  radial  forces  under  certain  modes  of  op¬ 
eration.  Now  the  currents  are  allowed  to  automatically  adjust 
themselves  according  to  the  relative  airgaps  around  the  ma¬ 
chine.  The  same  firing  angles  from  the  previous  run  were  used. 
The  radial  force  is  plotted  in  the  Fig.  11  and  is  observed  to  be 
reduced  from  that  of  the  current  mode  case.  An  average  radial 
force  of  88  lbs  is  obtained  from  the  MEC  dynamic  analysis.  It 
should  be  noted  that  even  if  the  external  terminals  were  ex¬ 
cited  with  a  current  source,  the  currents  would  still  divide 
asymmetrically  according  to  the  unbalance  in  inductances  if 
the  poles  are  connected  in  parallel.  An  natural  flux  balancing 
would  still  occur  and  the  radial  forces  would  be  reduced  from 
the  series  connected  case.  Generally  speaking,  poles  are  con¬ 
nected  in  parallel  rather  than  series  to  conform  to  the  voltage 
and  current  constraints  of  the  power  electronics  that  typically 
favor  lower  voltage. 

VII.  Radial  Force  Frequency  Response 

A.  Motivation 

So  far  the  dynamic  analyses  indicate  that  the  radial  force 
appears  as  a  negative  stiffness  in  the  direction  of  the  eccentric¬ 
ity.  However,  the  eccentricity  has  only  been  treated  statically. 
That  is  the  rotor  has  been  offset  from  its  centered  position  and 
then  spun  about  its  geometric  center.  When  integrated  with  a 

magnetic  bearing  system,  the  eccentricity  will  also  be  a  d^ain- 

ic  variable.  In  the  case  of  a  whirl  condition,  the  eccentricity  will 
be  a  rotating  vector.  The  whirl  may  be  forwards,  backwards, 
synchronous,  sub-synchronous  or  super  synchronous.  The 


6 


18/12  SR  13,000  RPM.  250kW.  33%  eccentricity 


Rotor  Mechanical  Position  (degrees) 

Fig.  1 1  Dynamic  radial  force  comparison  for  voltage  and  current 
source  excitation. 


turbations  in  flux  and  one  radial  direction,  hence  the  other  two 
variables  will  be  dropped  from  the  notation. 


(12) 


Linearization  of  the  preceding  equations  results  in  the  fol¬ 
lowing  small  signal  equations  that  govern  the  behavior  for 
small  excursions  about  an  operating  point. 


Av  =  a- 


i-AXR+  d- 


'-AxR+- 


(») 

fiJ^A 

Finally,  the  transfer  function  from  the  position  to  the  radial 
force  is  derived  as  follows: 


(15) 


whirl  can  also  be  a  harmonic  combination  of  various  frequen¬ 
cies.  The  realistic  expectation  of  these  scenarios  raises  the 
question  of  the  frequency  response  for  the  radial  force.  While 
the  radial  force  can  essentially  be  mapped  as  a  negative  stiff¬ 
ness  as  a  function  of  torque-speed  operating  point,  the  fre¬ 
quency  response  of  the  force  will  be  of  primary  mterest  frorn 
the  rotordynamics  point  of  view.  The  influence  of  whirl  an 
stability  will  be  the  key  motivations  regarding  rotordynamics 
analyses.  Considering  awhirl  condition,  the  ecantricity  of  the 
rotor  is  changing  with  time.  In  any  one  radial  direction 
centricity  will  be  a  harmonic  sinusoidal  function  in  time.  Thus 
it  is  important  to  consider  the  response  of  the  radial  for^  to 
perturbations  in  the  eccentricity.  The  capability 
analysis  again  indicates  the  versatility  and  power  of  the  MBC 
model  coupled  to  a  time  stepping  electrical  and  mechanical 
simulation. 

B.  Simple  Analytic  Model 

Without  proceeding  to  the  full  MEC  model,  a  good  indica¬ 
tion  of  the  response  can  be  derived  from  some  simple  physical 
equations  and  some  basic  assumptions.  Since  the  force  is  de¬ 
pendent  upon  flux  established  in  the  magnetic  arcuit,  the  fre- 
quenty  response  of  the  radial  force  is  expected  to  essentially 
be  an  RL  response.  A  circular  whirl  will  appear  as  two  ortho¬ 
gonal  (sin  &  cosine)  functions  in  the  x  &  y  radial  directions. 
Considering  one  eccentricity  direction  only,  neglecting  the 
leakage  flux  and  eddy  currents,  and  assuming  radial  flux  in  the 
overlap  region  only,  a  few  physical  equations  can  quic^  re¬ 
veal  the  form  of  the  radial  force  frequency  response.  These 
equations  are  derived  as  follows  for  a  single  magnetic  path  m 
the  electrical  machine: 


A  =  X{i,x,y,e)  and  i  =  i(A,x,y,0)  (H) 

Note  that  the  flux  linkage  is  in  general  a  function  of  four 
variables.  The  following  equations  will  only  consider  per¬ 


il  is  evident  from  (15)  that  the  radial  force  frequency  re¬ 
sponse  is  essentially  the  same  as  that  of  the  flux  linkage  or  cur¬ 
rent  of  an  RL  circuit,  where  is  the  incremental  indue- 

tance  reciprocal. 

"While  this  result  is  quite  simple,  it  should  be  kept  in  mind 
that  it  represents  a  linearized  soultion  about  an  operating 
point  for  just  two  pole  pairs.  Considering  the  whole  machine, 
there  will  me  a  much  more  complex  magnetic  circuit  active  in¬ 
cluding  multiple  pairs  of  excited  poles.  Additionally,  a  whirl 
condition  may  have  a  considerably  sized  amplitude  that  ex¬ 
ceeds  the  small  signal  approximation  of  this  analysis.  A  de¬ 
tailed  simulation  of  the  machine  in  such  a  condition  is  re¬ 
quired  to  rigorously  address  the  frequency  response 
properties  for  eccentricity  variations.  None— the  less,  reason¬ 
able  results  are  obtained  from  this  simple  model.  Additional¬ 
ly,  the  extremes  of  the  response  can  be  quite  accurately  as- 
sessed  as  given  in  the  next  section. 

C.  Frequency  Response  Predictions 

A  simulation  model  based  on  the  machine  parameters,  flux 
linkage  characteristics,  and  the  analysis  of  the  preceding  sec¬ 
tion  is  used  to  investigate  the  frequency  response  of  the  radial 
force.  Specifically  a  perturbation  is  given  in  the  eccentnaty 
and  the  radial  force  is  taken  as  the  output.  Several  operating 
points  on  the  flux  linkage  curve  are  investigated  to  reveal  the 
rough  envelope  of  response  characteristics.  Fig.  12  illustrates 
the  simulation  results  for  various  operating  points  ranging 
from  maximum  to  minimum  incremental  inductance.  The 
magnitude  and  phase  of  the  complex  dynamic  stiffness  are 
plotted  in  typical  Bode  format.  Additionally  as  common  to  ro- 
tordynamic  analyses,  the  real  and  imaginary  parts  are  plotted 
as  stiffness  and  damping.  Note  that  the  imaginary  part  is 
scaled  by  frequency  to  give  equivalent  damping.  It  should  be 
observed  that  the  inherently  high  inductance  and  low  resist¬ 
ance  of  the  machine  result  in  substantially  low  break  frequen¬ 
cies.  This  indicates  that  the  SR  radial  forces  do  not  act  purely 
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angle  (F/x)  (degrees)  IF/xJ  (per  unit) 


as  a  negative  stiffness  in  the  direction  of  the  rotor  offset  for  th 
case  of  dynamic  eccentricity,  such  as  whirl. 


Fig.  12  Radial  force  frequency  response  to  eccentricity  variations. 


D.  Rotordynamic  Implications 

The  unbalanced  radial  force  is  of  certain  importance  to  the 
rotordynamic  consideration  of  the  mechMical  rotor  system. 
The  radial  force  presents  additional  loading  to  the  magnetic 
bearings  that  must  be  sized  to  account  for  the  force  if  not  rom- 
pensated.  Additionally,  the  radial  force  is  a  negative  stiffcess 
which  is  a  destabilizing  property.  The  pulsating  nature  of  the 
radial  force  as  the  machine  is  operated  dynamically,  presents 
the  potential  for  excitation  of  various  resonant  modes.  Since 
the  frequency  is  synchronous  to  the  spin  speed  and  36  times 
higher,  it  varies  continuously  and  will  inevitably  pass  throu^a 
wide  range  of  frequencies  and  critical  speeds  of  the  rotor.  TOe 
frequency  response  as  discussed  in  the  previous  section  m- 
plies  a  phase  shift  in  the  radial  force  under  dynamic  eccentrici- 
n  conditions,  namely  whirl.  The  phase  shift  in  turn  implies  an 
effective  damping.  This  damping  has  the  potential  to  be  posi¬ 
tive  (stabilizing)  or  negative  (destabilizing).  The  RL  natme  of 
the  time  constant  suggests  that  the  response  will  be  typically 
very  fast  as  compared  to  the  time  constants  of  the  mechamcal 
system.  None-the-less,  the  complexities  associated  with  the 
highly  non-linear  behavior  of  the  magnetics,  machine  operat¬ 
ing  modes  (i.e.  commutation,  etc.),  and  complex  rotordyn^ic 
interactions  maintain  significant  challenge  in  understanding 
the  behavior  of  the  system  and  potentially  in  the  satisfactory 
Operation  of  the  system. 

VIII.  Radial  Force  and  the  Torque-Speed 
Operating  Range 

It  is  of  interest  to  map  out  the  radial  force  as  a  function  of 
torque  speed  operating  point  for  the  starter/generator.  Fig.  13 
illustrates  the  torque-speed  operating  envelope  of  the  start¬ 


er/generator.  Superimposed  on  the  plot  ar  reiaffications  of  the 
average  radial  force  for  various  points  of  (.operation. 


Fig.  13  Radial  force  -vs-  Torque-speed  operati.ingpomt,  18/12 
SRM,  33ft  eccentricity. 

IX.  Experimental  Radial  Force  MIeasurements 

To  gain  further  confidence  an3  validatiion  for  the  radial 
force  predictions,  a  commercially  available  lV2/8test  motor  was 
experimentally  investigated.  The  motor  is  (Considerably  small¬ 
er  than  the  18/12  machine  and  is  fabricated  ffrommore  conven¬ 
tional  materials.  Simple  pull-force  tests  wcerc performed  pn 
the  test  motor  to  generate  the  force-vs-currrent  characteris¬ 
tic  plots  as  predicted  for  the  18/12  machine, .  A  finite  element 
model  of  the  test  machine  was  developed  ba^  on  measured 
dimensions  and  on  an  estimated  magnetic  miatcrial.  Fig.  14  il¬ 
lustrates  the  experimental  and  analytical  ressults.  The  results 
confirm  the  qualitative  nature  of  the  force— ws-cutxent  char¬ 
acteristic.  The  discrepanc)’  in  absolute  accuirary  is  most  likely 
due  to  the  measurement  procedure  and  to.tincxact  BH  data. 
While  additional  experimental  results  are  n(bt  available  at  this 
time,  the  finite  element  results  have  been  useiful  in  the  verifica¬ 
tion  of  the  MEC  model. 


Phase  Current  (amp) 

Fig.  14  Predicted  and  measured  radial  force  -vs-  current  for  a  12/8 
SR  test  machine. 
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X.  Conclusions 

A  powerful  magnetic  equivalent  circuit  model  has  been  de¬ 
scribed  and  employed  to  investigate  the  static  and  dynmic  un¬ 
balanced  magnetic  radial  force  properties  of  a250kW  integral 
starter  generator  switched— reluctance  machine  operating 
with  asymmetrical  mechanical  rotor  positions. 

The  static  unbalanced  radial  forces  as  a  function  of  excita¬ 
tion,  rotor  angular  position,  and  eccentricity  has  been  pre¬ 
dicted  and  validated  with  finite  element  results.  The  magnetic 
model  was  then  integrated  with  temporal  differential  equa¬ 
tions  to  produce  a  complete  time  stepping  simulation  of  the 
SR  drive  system  with  the  capability  of  modelling  asymmetrical 
rotor  operation. 

The  dependency  of  radial  forces  on  the  torque-speed  op¬ 
erating  point  as  well  as  excitation  mode  and  winding  configu¬ 
ration  has  been  investigated.  Specifically,  it  has  been  deter¬ 
mined  that  voltage  mode  excitation  of  the  machine  with 
parallel  connected  windings  has  inherent  flux  balancing  and 
results  in  substantially  reduced  radial  forces  (by  a  factor  of 
three  or  four)  from  that  of  series  connected  windings.  Finally, 
limited  experimental  results  with  a  12/8  test  machine  confirm 
the  static  analytical  predictions. 

The  model  shows  that  the  expected  level  of  radial  forces  is 
such  that  it  will  be  necessary  to  utilize  some  means  of  radial 
force  compensation  in  a  magnetic  bearing  switched-reluc- 
tance  system.  Future  work  will  utilize  the  model  to  investigate 
the  performance  of  alternative  compensation  techniques. 
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